Possible complex annihilation and B Kn direct CP asymmetry 
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■ We point out that a sizable strong phase could be generated from the penguin annihilation in the 

' soft-coUinear effective theory for B meson decays. Keeping a small scale suppressed by 0{A/mt), A 

being a hadronic scale and mt the b quark mass, in the denominators of internal particle propagators 
without expansion, the resultant strong phase can accommodate the data of the A'^tt^ direct 

CP asymmetry. Our study reconciles the opposite conclusions on the real or complex penguin 
annihilation amplitude drawn in the soft-coUinear effective theory and in the perturbative QCD 
' approach based on kr factorization theorem. 
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Qh| The effect of scalar penguin annihilation on charmless nonleptonic B meson decays has attracted intensive 

pi • attention. This power-suppressed contribution is chirally enhanced, i.e., proportional to ^p/rrib in B ^ PP 
D \ decays, where jip is the chiral scale associated with the pseudoscalar meson P and the b quark mass. Since 
j-^^' it involves endpoint singularities, it was parameterized as a free parameter Xa = lii(mf,/A)[l + p^exp(i0yi)] 
in QCD-improved factorization (QCDF) [jy], with A being a hadronic scale, and pA and (j)A varied arbitrarily 
PsJ ■ within some artificially specified ranges. In order to fit data such as the i?° —> K^tt^ direct CP asymmetry 
^ I Acp(-B*' — > K^TT^), (f>A must take a sizable value. On the other hand, the contribution from scalar penguin 
annihilation has been found to be almost imaginary in the perturbative QCD (PQCD) approach based on kr 
factorization theorem [3, 0] , and the resultant strong phase leads to a prediction consistent with the measured 
^N^j ■ Acp(-B° — > K^-K^). The annihilation amplitude was not considered in the leading-power formalism of soft- 
collinear effective theory (SCET) . Instead, a nonperturbative complex charming penguin was introduced 

,^ ^ to accommodate the data of Acp(i3° A'^tt^). In the recent SCET formalism with the zero-bin subtraction 

\ 0' annihilation contribution becomes factorizable, and has been concluded to be almost real [s*]. 
. The motivation of this paper is to reconcile the opposite theoretical observations on the almost imaginary or 

almost real penguin annihilation derived in PQCD and in SCET. We shall first point out that the comparison of 
the measured Acp (i?^ — > -ftT^Tr^) and Acp(i?^ K^p^) indicates an imaginary penguin annihilation amplitude 
\ d The B^ -> K^n° {B^ K^P°) decays involve a. B ^ P {B ^ V) transition, so the penguin 
^ . emission amplitude is proportional to the constructive (destructive) combination of the Wilson coefficients 
a4 -|- (— )2(/ii^/TO;,)a6, pLK being the chiral scale associated with the kaon. The annihilation effect is then less 
inffuential in the former than in the latter. If the penguin annihilation is real, both decays will exhibit small 
direct CP asymmetries, i.e., Acp{B^ — + if^7r°) Aqp{B^ — > K^p'^) w 0. If imaginary, it will cause a larger 
Aip(B± ^ K^p°). The current data ^cp(S± K^tt^) = 0.050 ± 0.025 and Acp(B± ^ 0.3ll|J;io 
[111 ] favor an imaginary penguin annihilation. 

We emphasize that strong phases, generated by subleading corrections, are the leading effect for direct CP 
asymmetries of B meson decays. For example, the prediction for the direct CP asymmetry Ac-p{B^ if^Tr") 
is sensitive to the strong phase of the ratio C/T [ij, where C (T) is the color-suppressed (color- allowed) 
tree amplitude, though the branching ratio B{B^ K^tt^) is not. Assuming this ratio to be real as in the 
leading-power SCET it is difficult to explain the data. Therefore, the study of strong phases requires a 
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careful treatment of subleading corrections. It will be explained that the different penguin annihilation effects 
observed in PQCD and SCET arise from whether parton transverse momenta fcy and other intrinsic mass scales 
in particle propagators are expanded or not. If these small scales are neglected or expanded, the internal particles 
in an annihilation amplitude are on their mass shell only at the endpoints of parton momentum fractions, where 
hadron distribution amplitudes usually vanish, or the zero-bin subtraction applies. An annihilation amplitude is 
then real. Including kx, the on-shell condition of internal particles does not occur at the endpoints, so that there 
is a potential to generate a sizable strong phase. We claim that when to;, approaches infinity, the on-shell region 
coincides with the endpoints, and the same vanishing results for strong phases will be derived, irrespective of 
whether the small scales are expanded into a power series. For the physical value of m^, however, a formally 
power-suppressed correction may have a significant numerical effect on strong phases, and lead to large direct 
CP asymmetries in B meson decays. 

As argued in Ref. [l4[ , a parton, carrying a transverse momentum kx as small as a hadronic scale A initially, 
accumulates its kx after emitting infinitely many collinear gluons. When the parton participates in a hard 
scattering eventually, kx can become as large as the hard scale. Such an accumulation is described by the 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi evolution [l5| for a parton distribution function in inclusive processes 
and by the Sudakov evolution for a hadron wave function in exclusive processes. For two-body nonleptonic 
B meson decays, k^ of internal particles in a hard kernel reaches the hard scale of 0{mbA). That is, the effect 
resulting from kx is suppressed by a power of r = k'^/ml ^ 0{A/mb)- In SCET, the power counting rule for 
kx is different, which is always treated as being 0(A), and expanded. However, there exists a scale of O(TObA) 
from the hard-coUinear modes, which is also suppressed by A/to^ compared to to^. To verify the above claim, 
we shall keep a small scale in particle propagators, which can be regarded as an averaged parton transverse 
momentum in PQCD or the hard-coUinear scale in SCET, and examine its effect on the penguin annihilation 
in the SCET formalism with the zero-bin subtraction (3]. 

Before computing the direct CP asymmetry of the K^tt^ decays, we illustrate why a formally power- 

suppressed correction of 0(r) could produce a sizable strong phase in an annihilation amplitude. Expand a 
kernel of the form 

^T--4-+^(-)' 

X — r + l€ X + l€ \x/ 

which appears in a convolution with a meson distribution amplitude. Eq. ([1]) holds in principle as long as the 
contribution from the small x region is suppressed by the meson distribution amplitude, namely, as the main 
contribution comes from the region with r/x <^ 1. On the other hand, we have the principle- value prescription 
without expansion, 

^ = P— inS{x - r) . (2) 



X ~ r + ie 

Convoluting the kernel with the distribution amplitude = 6x(l — x), the real parts from Eqs. ((T|) and 
differ by only 15%. The imaginary part from Eq. ([T]) vanishes, but that from Eq. ^ reaches half of the real 
part for a typical value of r ~ A/to^ ~ 0.1. Obviously, in order that the imaginary part becomes negligible, i.e., 
about 5% of the real part, r must decrease to 0.01 (or nib increases up to 50 GeV). The lessons we learn from 
this simple example are 1) as x has the substantial probability to be close to r, which is small but away from 
the endpoint, the expansion in a power series of r breaks down, and an imaginary piece could develop; 2) the 
expansion is reliable only for sufficiently small r such that the contribution from x ^ r is highly-suppressed like 
the endpoint one; 3) r is expected to give a minor (larger) effect on branching ratios (direct CP asymmetries) 
of B meson decays. 

Let the momenta of the outgoing quark u and antiquark u in opposite directions be k2 = {0,yP^ ,0x) and 
ks = {xP^ ,0,0x), respectively, for the decay B^ K~n~^, where P2 (P3) is the pion (kaon) momentum and 
X — 1 — X. We quote the expression for the penguin annihilation amplitude in the SCET formalism with the 
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operator in the effective weak Hamiltonian. 



zero-bin subtraction 
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where Gf is the Fermi constant, fB,K,-K the meson decay constants, Xu.c the products of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements, fit ~ rub the hard scale, Ci the Wilson coefficients, and //jr the chiral scale 
associated with the pion. The logarithmic terms ln/i± in Ref. [8] have been dropped since they are cancelled 
by the corresponding logarithms in the convolutions. Because of the large theoretical uncertainty shown below, 
the constant k resulting from the above logarithmic cancellation will be neglected Q . The three-parton twist-3 
contribution to the penguin annihilation, being numerically smaller by one order of magnitude than Eq. ^ 
17| . is not included. 

Motivated by the illustration based on Eqs. ([T]) and 0, we introduce a small constant r into internal quark 
propagators involved in the factorizable piece of Eq. Q, corresponding to Fig. [T] Inserting r into internal 
gluon propagators generates a strong phase down by a factor three. The strong phase from the nonfactorizable 
annihilation amplitude is smaller by two orders of magnitude. We stress that adding r in the aforementioned 
way causes a double counting of the contributions from higher-order operators in SCET, and should be regarded 
as only a test of our claim. Applying the principle-value prescription, we obtain the extra imaginary pieces via 
the following substitutions, 
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Employing the parameterizations for the leading-twist distribution amplitude (puix) and for the two-parton 
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twist-3 distribution amplitudes (j^^ix) [1, 17 1 

(j)Mix) = 6x{l - a;) [l + af (6x - 3) + 6af (1 - 5x + Sx^) - lOaf (1 - 9x + 21a;2 - Ux^) 

+15af (1 - 14a; + 56x'^ - 84x^ + 42a;'') H ] , 

0^(a;) = 6a;(l - a;) [l + at^p(6x - 3) + 6a^^p(l - 5x + Sx^) + •••], (6) 

with M — TT, ii', it is easy to find that both Eqs. ^ and (O are proportional to r as expected. 

The importance of the penguin annihilation contribution relative to the full penguin one has been estimated 
in SCET [8(, and found to be about 10% with large uncertainty in the K^tt^ decays. The full penguin 

contribution does not come from an explicit evaluation in the same SCET framework, but from a fitting to the 
B Kn data. We can certainly follow this approach. However, the factorization formulas for the emission 
amplitudes have been available in Ref. Q , so they will be adopted in the numerical analysis below. The feature 
of generating strong phases does not depend on how we estimate the emission amplitudes. Besides, we shall 
not include the free parameters associated with the long-distance charming penguin, which is not factorizable 
in SCET. As demonstrated later, a decay am plit ude under the zero-bin subtraction is very sensitive to higher 



Gegenbauer moments a^^ and in Eq. ([6]) [18|, which are mostly unknown. Hence, we shall determine these 
moments by fitting the SCET formulas to data of branching ratios, which are then used to predict direct CP 
asymmetries. If a strong phase from the source considered here is sizable, the whole CP asymmetry cannot be 
attributed to the nonperturbative charming penguin alone. 

At lowest order in asimt) with the Wilson coefficients r(+' = 1 and C^+^ 1 in SCETi fi|, the S ^ tt 
transition form factor is decomposed into 

U{E)^C''''{E)+Cf^E). (7) 

The second term is factorizable, written as 



. . Mrns 4^^ . 2^ 2E . ^ ^ 

^ ' AE-^ 9 Vms rrik " V Jo k+ 



(8) 



where ~ V'ti&A is the intermediate scale, and fc+ the momentum of the spectator quark in the B meson. For 
charmless two-body nonleptonic B meson decays, we take the pion energy E = ttib/'^, ms being the B meson 
mass. The first term also becomes factorizable after implementing the zero-bin subtraction for the endpoint 
singularity 0], 

.B.r^^_ fB^niB ^naM f\^. r f i^ + y)My) , .. (0^ + ^0(2/) 



where only the terms from the two-parton pion distribution amplitudes are kept. The relation among 0^ 
and can be found in Ref. The formulas for the B ^ K form factor in SCET arc similar. We multiply 
Eq. ([7j) by the appropriate CKM matrix elements and Wilson coefficients, including a part of ncxt-to-leading- 
order corrections [19], to obtain the emission contributions from both the tree and penguin operators. The 
Wilson coefficient ag was neglected in the previous SCET analysis, since the associated penguin contribution 
is power-suppressed. However, it is enhanced by the chiral scale, and numerically crucial. Furthermore, the 
power-suppressed annihilation has been formulated into SCET, so there is no reason for ignoring ae [l9l |. 

The zero-bin subtraction for the logarithmic endpoint singularity associated with the pion distribution am- 
plitude (t>Tr in the ffi-st term of Eq. ([9]) is referred to Ref. [7], where the term proportional to y in (1 -|- y) does 
not require subtraction. We also need the zero-bin subtraction for the linear endpoint singularity present in the 
second term of Eq. ([91) 



2£ 
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where n ■ P2 = 2E. The subtraction associated with the derivative of the two-parton twist-3 pion distribution 
amphtude, (j)'^' , is similar. 

We consider the models for the B meson distribution amplitudes (f)g proposed by Kodaira et al. (KKQT) 



2ll | and by Grozin and Neubert (GN) [22|. The associated zero-bin subtraction is defined by 



> . Hk- _ /V + m r^Msio) , (11) 



(fc+)0 Jq k+ Jo k+ V ^ 

-^(1 -ln2) + ln(— )0b(O), for KKQT 



A _ 

if , A \ , / n • u A N 



- 7i5+ln- +ln(^^)05(O), forGN 



^ r*-4§r^- +lnf^)4'(0), (12) 



(fc+2)0 - fc+2 A;+ ^ V ^+ 

_i_ In 2 + In (!^) 0+'(O) , for KKQT 

-^f7£; + ln-i+ln(^)0+'(O), forGN 

with the parameter relation luq = 2A/3, A being the B meson and b quark mass difference. In the above 
expressions n is a light-like vector along the Wilson line in the definition for the B meson distribution amplitudes, 
and V is the B meson velocity. The terms containing ln/i± in Eqs. p0p -()12 |l are also dropped. 

For the numerical analysis, we assume the Gegenbauer moments of the pion and kaon distribution amplitudes, 
= 0.0, af = -0.05 consistent with the results in Ref. [H,!!!!, af = = 0.2 [Ulli,!!!!, = 0, af = aj. 



^ipp ~ "^ipp ~ ^'^'^ '^^pp ~ "^^pp! among which a*^ and are most uncertain. To simplify the formulas, we 
do not consider the Gegenbauer moment af for the twist-2 kaon distribution amplitude. That is, we keep one 
most uncertain parameter from each of ipM and 0^^, whose variation is sufficient for our purpose. The hard and 
intermediate scales are fixed at /i^ — nib and fii = i/Amf,, respectively, with A = 0.55 GeV and nib = "m-l^ = 4.7 
GeV. Other relevant heavy-quark masses are taken to be nic = = 1.4 GeV and nib = rn^^{rnb) = 4.2 GeV. 
We obtain the chiral scales fiTrifih) = 2.4 GeV, ^Kip-h) — 3.0 GeV, fJ.-K{fJ,i) = 1.8 GeV, and ^MKifJ-i) = 2.3 GeV 
from the two-loop running for the strong coupling constant with as{Mz — 91.1876 GeV) = 0.118 and for the 
light-quark masses with niu.d{'2 GeV) = 5 MeV and ms{2 GeV) = 95 MeV. We take the Wilson coefficients for 
four-fermion operators evaluated at fih = m-b and at next-to-leading-logarithmic level: Ci = 1.078, C2 — —0.177, 
C3 = 0.014, C4 = -0.034, C5 = 0.009, Cq = -0.040, C7 = 0.7 x 10-^ Cg = 4.5 x lO"*, Cg = -9.9 x lO'^, 
and Cio = 1.8 x lO"'^. Those for dipole operators at leading-logarithmic level are Cr^ = —0.314 and Csg = 
-0.149 [27|. We also take the Fermi constant Gp = 1.16639 x 10"^ GeV'^, the decay constants fs = 0.22 
GeV, Jk = 0.16 GeV, and U = 0.131 GeV, the meson masses ms = 5.28 GeV, niK = 0.497 GeV, and 
= 0.14 GeV, the B meson lifetime = 1.530 x 10~^^ sec, and the CKM matrix elements Vus = 0.2257, 



Vub = (4.2 X IQ-^) exp(-i(/)3), Vcs = 0.957, and Vcb = 0.0416 with the weak phase ^3 = 74° [28| 



Adopting the above parameters, the two pieces ( and C,/ of the B ^ tt form factor are written as 

^ r 0.01 + 0.75 aj + 2.57 + 0.43 ajpp, for KKQT 
^ ^ \ 0.09 + 0.65 aj + 2.23 aj- 2.73 ajpp, forGN 

^ r 0.016(1.0 + aj + aj), for KKQT 
^■^ \ 0.024(1.0 + aj + aj), forGN. ^ > 

Note that the coefficients in Eq. (fT5|) grow quadratically with the order n of the Gegenbauer moments [3] ■ 
This sensitivity is attributed to the increasing slope of the higher Gegenbauer polynomials at the endpoints 
of the momentum fraction x. The sign flip of the ajpp terms indicates that C^^'^ also depends strongly on the 
models of the B meson distribution amplitudes in SCET. We mention that the PQCD approach does not suffer 
such sensitivity, because the endpoint singularity is smeared by including parton transverse momenta /ct, whose 
order of magnitude is governed by the Sudakov factor. 
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The strong dependence on the higher Gegenbauer moments also appears in the penguin annihilation ampli- 
tude, 



\/2 ALanniK-7r+) 

K 



lo^pr" = -10 



FTO^ (IGeV) 

= 2.76(0.07 + aj ) (1.20 + aj ) + a2pp(27.0 + 413.1 aj) -iirr ajpp(53.2 + 1747 aj) , (15) 

with a significant growth of the coefficients of aj. The imaginary contribution is proportional to the second 
moment flj^p. In fact, it could depend on the zeroth moment, i.e., the normalization of (p^, if the denominator 
1 — y is not replaced by 1 in the subtraction term in Eq. The denominators 1 — y and 1 correspond to 
different zero-bin subtraction schemes. 

Note that the size of the imaginary part depends on the amount of the subtracted contribution, i.e., on 
zero-bin subtraction schemes, since it is generated at i ~ A/mf, ot y ^ A/m^ as shown in Eqs. ^ and ([5]). 
The dependence on subtraction schemes also exists in all other definitions like Eqs. (fTUl) - P^ . which will not be 
discussed in this work. 

For the range of aj, the crude bound aj > —0.07 has been determined in Ref. |2fl]. The analysis based on the 
data of the pion transition form factor suggests aj « —0.05 in Ref. [29] and the constraint -l-aj — — 0.03±0.14 
in Ref. [s^l, both of which prefer a negative value of (considering a J « 0.2). The range of ajp^ is basically 
undetermined. We shall regard these two parameters as being free, and fix them by the strategy stated before: 
Adjust aj and aj^p, such that the B n form factor has the value around /+ = 0.24 ± 0.05 [lH, and the 
B" ^ K^t:^ decays have the branching ratio close to the data B{B° -> K^i:^) = (19.4 ± 0.6) x 10"^ [ll]. 
Because the last two terms in Q^'^ for the KKQT model are of the same sign, and the coefficient of a\ is large, 
the constraint from the form factor value leads to a smaller aj. Eq. (fT5|) then implies that the coefficient of r, 
i.e., the imaginary part of the annihilation amplitude, is smaller, and that the strong phase is less sensitive to 
the variation of r. On the contrary, the last two terms in Q^'^ for the GN model have the coefficients with the 
same order of magnitude, but in opposite signs. Hence, aj (and also ajpp) is larger, and the strong phase is 
more sensitive to the variation of r in this case. 

Employing the KKQT model for the B meson distribution amplitudes, we obtain a\ m 0.01 and Ojpp ~ 0.23, 
corresponding to which the B ^ n form factor, the B° K^n^ branching ratio, and the predicted direct CP 
asymmetry are given by 
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B{B° 



AcpiB° K'^n^) = { 0.05 for r = 0.1 (16) 



We do not attempt a fine tuning here, but accept the values of aj and ajpp as solutions, when they produce the 
B ^ TT form factor and the 5° K^n^ branching ratio close to the designated ranges. The results shift with 
the slight variation of aj and ajpp, but the behavior for different r in Eq. (jl6p has the same pattern. In principle, 
(^'^ and (j'^ have the same scaling law in ag and in 1 /m^ f3l], 'sij . The numerical hierarchy (^'^ ^ (j'^ in 
Eq. (fT6|) . consistent with the PQCD results may be altered in different zero-bin subtraction schemes. It is 
obvious that the power correction associated with r has a negligible effect on the branching ratio. However, the 
power correction generates a strong phase: Acp{B^ K^n^) decreases by 40% from r = to r = 0.1. Since the 
imaginary part is proportional to r, it is difficult to accommodate the data Acp{B^ — > K^n^) — — 0.097±0.012 



llj l with a reasonable value of the power-suppressed r using the KKQT model. 
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For the GN model, we find two sets of solutions corresponding to aj 0.18 and a2„„ ~ 0.15 
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2pp 

= 0.21, 
(f"" = 0.03, 

r 20.1 X 10-6 for r = 0.0 
B{B° -> K^TT^) = < 20.4 X 10-6 for r = 0.1 
[ 25.1 X 10-6 for r = 0.2, 

( 0.06 for r = 0.0 

Acp{B" K^TT^) = < -0.06 for r 0.1 (17) 
[ -0.14 for r = 0.2, 

and to aj w -0.22 and aj^p w -0.20, 

C^^ = 0.28, 
Cf" - 0.02, 

Acp(S° ^ i^T^TT^^) I -0.10 for r = 0.1 (18) 
[ -0.20 for r = 0.2 . 

The existence of the two sets of solutions with opposite signs is understandable. Because the term proportional 
to aj in the imaginary part of Eq. dominates over the constant term as |aj| reaches about 0.2, the product 
aJppOj matters, and aj and ajpp can flip sign simultaneously. 

As indicated by Eqs. (fT7|) and (fT8|l . the branching ratio is stable, while the strong phase is very sensitive 
to the variation of r, so that we easily accommodate the data of Acp{B'^ K^ir^) with a typical value of 
r = 0.1 ~ 0.15. The predicted Acp{B'^ — > K^n^) for r — 0, i.e., real penguin annihilation (r = 0.1, i.e., complex 
penguin annihilation) is close to that from QCDF in the default scenario [33] (PQCD 0,111). Therefore, the 
strong phases resulting from the power-suppressed source in the penguin annihilation could be numerically 
crucial for the estimation of direct CP asymmetries. We then understand the opposite conclusions on the 
effect of the penguin annihilation drawn in SCET and in PQCD: The almost real annihilation amplitude in the 
former and the almost imaginary annihilation amplitude in the latter are attributed to the different treatments 
of the formally power-suppressed terms at the physical b quark mass. Note that the solutions of aj and ajpp in 
Eqs. (fT6|) -(fT8 | will be changed, if higher Gegenbauer moments in Eq. ([6|) are taken into account, which cause 
even larger variation of the decay amplitudes. However, the strong dependence of AcpiB'^ K^tt^) on r will 
persist. 

SCET provides a systematical expansion in powers of A/m?,, which is somewhat twisted here by keeping 
subleading terms in particle propagators in order to demonstrate a possible mechanism for generating strong 
phases. This twist of SCET actually violates its power counting rules and other aspects. Hence, our analysis 
does not imply the breakdown of SCET in its application to B meson decays, but helps clarifying why there 
are discrepancies in the study of direct CP asymmetries from SCET and PQCD. It hints that more caution is 
necessary for fixed-power evaluations of direct CP asymmetries at the physical mass mt. The expansion would 
be reliable for decay rates and direct CP asymmetries, if the b quark mass was 10 times heavier. In that case, 
the contribution from the on-shell region of internal particles can be really suppressed by hadron distribution 
amplitudes, or excluded by the zero-bin subtraction. For nib ~ 5 GeV, a novel method might be demanded. 

We have shown that introducing a small scale into denominators of internal quark propagators can accommo- 
date both the measured branching ratio and direct CP asymmetry of the B^ K^n^ decays. Keeping a small 
quantity in denominators without expansion is equivalent to resummation of the associated corrections to all 
powers. It is similar to resummation of part of higher-order corrections in for many QCD processes. It has 
been explained that at least the parton transverse momenta can be maintained in denominators consistently 
in kr factorization theorem [s^, This treatment is justified by different power counting rules, which hold 
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in the region of small parton momenta l35|. This alternative power expansion, postulated in kx factorization 
theorem, has led to strong phases in more agreement with the indication of data in B meson decays. 
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